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The  paper  deals  with  the  high-fidelity  modeling  and  simulation  of  a  powered  parafoil-payload  system  with 
respect  to  its  application  in  autonomous  precision  airborne  cargo  delivery.  In  the  proposed  concept  the  cargo 
transfer  is  accomplished  in  two  phases:  Initial  towing  phase  when  the  glider  follows  the  towing  vessel  in  a  passive 
lift  mode  and  the  autonomous  gliding  phase  when  the  system  is  guided  to  the  desired  point.  During  the  towing 
phase,  the  system  gains  as  much  altitude  as  possible  by  taking  the  angle-of-attack  that  will  provide  the  best  lift. 
Once  sufficient  altitude  is  attained,  the  gliding  phase  starts.  The  system  is  steered  to  the  desired  location  by 
controlling  the  lengths  of  the  rear  suspension  lines  using  two  control  inputs.  The  paper  presents  the  concept  of  the 
system,  its  6DoF  model,  the  control  algorithm  at  the  stage  of  passive  glide  and  the  simulation  results. 


I.  Introduction 

IN  the  past,  the  modeling  of  the  parafoil-payload  system  has  been  done  with  respect  to  various  applications  such  as 
NASA  assured  crew  return  vehicle,  payload  recovery  in  sounding  rocket  flights,  supply  troops  with  ammunition  and 
other  essentials  from  airplane,  etc.  In  all  these  applications  the  payload  along  with  the  parafoil  is  released  from  a  certain 
height  and  then  it  is  guided  to  the  pre-defmed  target  point  on  the  earth.  During  littoral  operations  there  is  a  demand  to 
supply  troops.  It  is  required  that  the  system  is  capable  of  delivering  goods  from  the  ship  onto  unprepared  spots  on  the 
shore,  be  passive  and  quiet. 

Figure  1  represents  the  underlying  concept  of  the  proposed  system.  The  550  sq  ft  double-skin  parafoil  with  750  lb 
payload  is  being  towed  and  released  at  the  altitude  of  about  2000  -  2500  feet.  It  then  gains  as  much  altitude  as  possible  in 
a  powered  mode  and  then  glides  shoreward.  It  is  assumed  that  the  airborne  autopilot  will  use  way-point  navigation  to 
deliver  supplies  to  the  predetermined  position  within  the  range  of  a  couple  of  miles.  Close  to  the  target  it  performs  a 
wind-determining  maneuver  and  safely  lands  into  the  wind. 

This  paper  describes  current  status  of  the  project,  developed  parafoil-payload  system  model,  mathematical  model  of 
the  system  and  results  of  simulations.  Along  with  the  analysis  of  passive  glide,  the  paper  focuses  on  the  influence  of  the 
powered  flight  on  the  dynamics  of  the  system.  The  paper  also  addresses  the  environmental  effects  on  the  dynamics  of  the 
flight  such  as  changes  in  the  wind  velocity  and  the  density  of  air.  It  finally  summarizes  the  observations  and  conclusions 
from  the  typical  simulation  results. 
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Figure  1.  Stages  of  deployment 


II.  The  Description  of  the  System 

Figure  2  shows  a  schematic  of  the  parasail-payload  system  labeled  with  the  4  major  points  used  in  the  equation  of 
motion  development.  The  tow  cable  emanates  from  a  location  near  the  stem  of  the  ship  at  point  1  and  is  attached  to  the 
parasail  system  at  point  3.  The 
tow  line  length  L{t) ,  and  its 

first  two  derivatives,  L{t) 

and  L{t)  ,  are  assumed  to  be 
known.  The  payload  and  the 
parasail  are  assumed  to  consist 
of  a  single,  6  degree  of 
freedom  rigid  body  with  mass 
m  and  inertia  matrix  I  .  Its 
center  of  mass  is  at  point  2. 

The  center  of  pressure,  where 
the  lift  and  drag  forces  are 
assumed  to  act,  is  at  point  4. 

Point  5  denotes  the  location  of 
the  constant  thrust  engine.  The 
inertial  frame  is  also  depicted 
in  Figure  1  and  will  be 
abbreviated  using  the 
notation {/}.  The  vector  from  developing  the  dynamic  equations  are  labeled. 

the  origin  of  {/}  to  the  ship’s  tow  point  is  denoted  as  and  is  assumed  to  be  known.  Thus,  the  tow  point’s  absolute 

velocity  and  acceleration  vectors  are  also  known  and  are  denoted  as  v^{t)  and  a^{t)  respectively. 
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III.  Modeling  of  the  System 

Figure  3  shows  a  close-up  of  the  parasail/payload  assembly  to  illustrate  the  body  frame,  used  in  the  dynamic 
equation  development.  Its  origin  is  at  the  center  of  mass  of  the  parasail/payload  system  and  is  oriented  relative  to  {/} 
using  a  standard  yaw,  pitch,  roll,  Euler  angle  transformation.  Thus,  the  rotation  matrix  that  transforms  a  vector 
represented  in  {/}  to  a  representation  in  \b},  denoted  as  /R  is 


R 


1  0  0 
0  cos^  sin^ 
0  -sin^  cos^ 


COS  9 

0 

-sin^ 

0 

1 

0 

sin0 

0 

cos^ 

cosy/  sin^^  0 
-sin^^  cosy/  0 


0 


0 


1 


(1) 


where  (j),  6  and  y/  are  the  roll,  pitch,  and  yaw  of  the  body. 


4  h 


The  translational  dynamic  equations  are  obtained  by 
applying  Newton’s  second  law  to  the  system,  considering 

aerodynamic  {F^),  gravitational  {F^),  tow  line  tension  forces 

{F^).  The  system  free  body  diagram  is  shown  in  Figure  4,  and 
the  resulting  translational  dynamic  equation  in  Eq.  2. 


■■f,  +  f^  +  f,  +  f^ 


(2) 


where  03  the  absolute  acceleration  of  the  system’s  center  of 
mass.  The  gravitational  force  vector,  represented  in  {/}  and 

denoted  as  ^F  is 


F  =-mgK. 


(3) 


Figure  3.  Close-up,  side-view  of  the  parasail/payload 
system  showing  the  body  frame,  angle  of  attack  ( CX), 
and  the  absolute  velocity  of  the  center  of  mass. 

The  force  due  to  the  constant  thrust  engine  acts  along  the 
direction  and  is  ^F^  =  FJj^  when  expressed  in{Zz}. 

The  aerodynamic  force  vector  is  further  decomposed  into 
lift  and  drag  components  as  shown  in  Eq.  3. 


Fa  —  Fiift  + 


The  magnitude  of  the  lift  force  is 


F 


lift 


(4) 


ClQlA  where 


is  the  lift  aerodynamic  coefficient,  is  the  lifting 
dynamic  pressure  and  A  is  the  effective  area  of  the 
parasail.  Similarly,  the  magnitude  of  the  drag  force  is 


Figure  4.  Free  body  diagram  of  the  parasail/payload 
system  illustrating  the  4  external  force  vectors. 


f 

drag 


■  C^QA  where 


Cjj  is  the  aerodynamic  drag 


coefficient  and  Q  is  the  total  dynamic  pressure.  Since  lift 
and  drag  are  caused  by  the  velocity  of  the  parasail  relative 
to  the  local  wind  velocity,  the  wind  relative  velocity  of 
the  center  of  pressure  point  is  introduced  and  defined  as 
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(5) 


at 

where  is  the  absolute  wind  velocity  in  the  vicinity  of  the  parasail  and  is  the  absolute  angular  velocity  of  as 
shown  in  Eq.  5  where  it  is  also  represented  in  {^}. 


= 

^cos^+(^cos^sin^ 

=  Q< 

0 

-0  +  ipcos9cos(l) 

¥ 

(6) 


where 

1  0 

-sin0 

Q  = 

0  cos^ 

cos^sin^ 

0  -1 

COS0COS(^ 

Since  it  is  assumed  that  the  parasail  is  fully  inflated,  the  location  of  the  center  of  pressure  relative  to  the  center  of  mass  is 
constant.  Therefore,  V4  can  be  simplified  as 

^4,.=^2  +  ^bX^42-K  (V 

The  lifting  dynamic  pressure  is  assumed  to  be  unaffected  by  the  span-wise  contribution  of  wind  speed  and  is  given  in  Eq. 

8. 


1 


Ql  ~  2^^  ^4lw,x'^  ^4/w,z 


(8) 


where  p  is  the  local  density  of  air  and  the  left  superscript  indicates  components  of  the  vector  in  \b},  specifically 


^4  /  w 


'^4lw,x 


'^4lw,y 


4/w,z 


(9) 


The  total  dynamic  pressure,  used  to  compute  is  simply 
6  =  ^/^4/w-V4/w 


(10) 


Both  the  lift  and  drag  aerodynamic  coefficients  are  assumed  to  be  linear  functions  of  the  angle  of  attack  and  are  given  in 
Eq.  11. 
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(11) 


~  ^Lfi  ^L,a^ 

~  ^Dfi  ^D,a^L 

where  q,  q,  and  „  are  assumed  to  be  constant  and  the  angle  of  attack,  (2,  is  the  angle  between  and 

the  wind  relative  velocity  vector,  as  given  by  Eq.  12. 


a  =  tan 


-  V 


4/v 


V  ^4/w,x  ) 


(12) 


The  direction  of  the  lift  vector  is  orthogonal  to  V^/n, ,  and  when  represented  in  {Zz},  is  given  by  Eq.  13. 

=  |F;,j^|[sinct  0  cosct/  (13) 


The  direction  of  drag  is  opposite  to  and  so  the  drag  force  can  be  represented  as  shown  in  Eq.  14. 

^dras  ~  1^4  /  iv  (14) 


where  is  the  unit  vector  associated  with  V4/„. 

Finally,  the  tow  line  force  is  approximated  as  being  along  the  line  connecting  the  ship  and  tow  point  on  the  parasail 
system  and  is  described  by  Eq.  15. 

Ft  =  \Ft%  (15) 


where  the  unit  vector  Tjj  is  a  function  of  the  ship’s  position,  which  is  known,  and  the  position  of  the  tow  point,  which  is 
a  function  of  the  parasail’s  states.  This  expression  is  given  by  Eq.  16. 


As  = 


A-A| 


(16) 


The  value  of  Fj.  is  an  unknown  but  will  be  solved  at  each  integration  step  when  the  accelerations  are  computed.  The 
details  of  this  will  be  described  shortly. 


The  rotational  dynamic  equations  are  derived  from  Euler’s  equation  and  the  free  body  diagram  of  Figure  3  and  are  given 
inEq.  17. 


Ia)+a)xFa)=  +  p^2^F^  +  p^2^F^ 


(17) 


The  terms  of  Eq.  17  have  all  been  introduced  above,  and  thus  is  completely  described. 
The  12  states,  integrated  during  the  simulation,  are  given  in  Eq.  18. 
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<t> 

e 

¥ 


(18) 


At  each  time  step  X  must  be  computed.  The  time  derivative  of  roll,  pitch,  and  yaw  is  easily  computed  using  Q, 
introduced  earlier,  and  extracting  from  the  state  vector,  specifically. 


<  ^  > 

¥ 


(19) 


The  next  element  of  X ,  ^V2,  is  readily  computed  by  extracting  from  x  as  shown  in  Eq.  20. 


(20) 


The  angular  acceleration,  and  the  time  derivative  of  solved  simultaneously  with  the  tow  cable  force 

magnitude,  Fj..  A  set  of  7  equations  in  7  unknowns  (three  components  of  three  components  of  d(^V2^/dt  and 

the  scalar  are  created  in  matrix  form  at  each  integration  time  step,  by  inverting  the  7x7  matrix,  the  terms  in  Eq.  1 8  are 
computed,  as  well  as  the  tow  cable  tension.  The  tow  cable  tension  is  not  used  in  the  ensuing  integration,  but  is  a  useful 
piece  of  information  from  an  analysis  perspective. 

The  first  3  equations  are  created  from  Eq.  1 7  where  they  are  written  below  terms  containing  solved-for  quantities  on  the 
left  and  known  terms  on  right. 


I&- ^32 =  -q}xF q)+  p,2  xF^+  p,2xF^ 


(21) 


The  next  3  equations  are  created  using  Eq.  2  and  noting  that 

d  V rj  If  ^  If ^ 

«2  =  -r^  ^2 

dt 

Thus,  the  quantity  required  for  the  time  derivative  of  Eq.  18  is 
m-^-F^p,2  =  -m  co,x  V2+  F^+  F^+  F^ 


(22) 


(23) 


The  last  of  the  7  equations  is  formed  by  taking  two  time  derivatives  of  the  tow  line  length  constraint  equation  of  Eq.  24 
Ai  •  Fb  =  A  (24) 
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which  gives 


V31  ■  ^31  =  LL{t)  (25) 

and 

^31  '  ^31  +  ^31  ■  ^31  -  LL  +  (26) 


Representing  all  quantities  in  \b},  and  grouping  unknowns  on  the  left  and  known  terms  on  the  left  converts  the 
acceleration  constraint  of  Eq.  26  to  the  last  equation  of  the  7x7  system  and  is  shown  in  Eq.  27. 


at 


■V3 


3=[- 


*«^x*V2+*aj- 


)1' 


PiX 


^^31  *^31  ^ 


(27) 


where 


d{%) 

(28) 


dCpx) 

dt 


+'’4x^3 


1 


It  should  be  noted  that  after  the  tow  line  is  released,  Eq.  27  is  omitted  from  the  acceleration  solution  procedure,  and  a 
6x6  system  is  solved.  Finally,  pitching  and  rolling  moments  and  increased  drag,  caused  by  asserting  the  parasail  brake 
lines,  are  applied  directly  to  the  simulation  as  external  moments  and  forces. 


IV.  Simulation  Results 


The  simulation  was  used  to  obtain  a  better  understanding  of  the  challenges  associated  with  tow-up  to  2000  feet.  The 
parasail/payload  system  needs  approximately  8  to  10  minutes  to  be  deployed  to  600  meters  (1969  feet).  With  a  towing 
vessel  speed  of  13  m/s  (25.3  knots),  the  towing  vessel  needs  to  cover  6240  to  7800  meters  (2.9  to  3.6  miles).  Wind 
conditions  may  modify  these  requirements.  Shorter  deployment  times  translate  into  larger  towing  forces,  which  are  one 
of  the  primary  limiting  factors  of 

deployment  time.  Wind 


To  illustrate  the  use  of  the 
dynamic  simulation,  a  tow-up  to 
landing  evolution  is  considered 
below.  The  sail  area  was  selected 
as  46.4  (500  ft^  )  and  the 

payload  weight  was  286  kg  (630 
lb).  The  cross  wind  was  set  at  2 
meters  per  second  (4.5  miles  per 
hour)  at  an  inertial  heading  angle  of 
30  degrees.  The  parasail  is  towed  in 
the  inertial  x  direction  until  it 
reaches  600  meters  with  cable 
being  let  out  at  3  meters/second 


Figure  5.  Simulated  maneuver  evolution  from  tow-up  to  landing. 
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(9.8  feet  per  second)  and  the  towing  ship  traveling  at  13  meters  per  second  (25  knots).  At  500  seconds,  the  tow  line  is 
released.  The  engine  is  turned  on  at  550  seconds  with  400  Newtons  of  thrust  in  the  body  x  direction.  At  600  seconds  the 
parasail  turns  into  the  wind  (a  30  degree  turn)  with  engine  cut-off  at  800  seconds.  The  landing  occurs  at  approximately 
960  seconds  into  the  flight.  The  landing  control  strategy  is  to  flare  up  the  pitch  angle  when  the  parasail  is  4  meters  above 
the  ground.  This  reduces  the  forward  motion  by  about  2  meters/second,  with  only  a  slight  gain  in  downward  motion 
(approximately  .3  meters/second).  During  the  flare,  the  pitch  angle  goes  to  approximately  20  degrees.  Larger  pitch  angles 
are  possible,  but  result  in  a  greater  increase  in  downward  speed. 

Figure  5  shows  the  entire  evolution  in  an  x,y,z  inertial  coordinate  system.  The  black  line  is  the  parasail  trajectory 
and  the  blue  line  is  the  towing  ship.  Figure  6  shows  the  line  tension  time  history,  which  abruptly  ends  at  the  release 
time.  This  is  kept  below  1600  N  using  the  parasail  control  system.  It  is  possible  that  a  combination  of  parasail  control 
and  a  tension  control  system  on  the  winch  may  yield  better  performance  from  a  time-to-release  perspective  and 
robustness.  Figure  7  shows  the  x-y  position  of  the  parasail  during  the  flight  and  Figure  8  shows  the  altitude  time  history. 
These  illustrate  the  change  in  heading  during  the  evolution  and,  of  course,  the  change  in  z  direction. 


Figure  6.  Tow  cable  tension  time  history  for  the  maneuver  considered  in  Fig.  5. 
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V.  Conclusions 

The  complete  6DoF  model  of  the  parafoil  system  has  been  developed  and  realized  within  the  MathWorks 
MATLAB/Simulink  environment.  The  model  demonstrates  an  adequate  response  for  the  control  inputs  and  matches  the 
result  of  the  flight  tests  well. 

An  active  control  strategy  is  crucial  to  successful  tow-up.  Two  methods  were  examined:  constant  tension  with 
variable  tow-line  feed  rate  and  attitude/brake  control  using  constant  tow-line  feed  rate.  Although  both  methods  appear  to 
work  in  simulation,  it  is  likely  that  some  combination  of  the  two  approaches  will  be  best  in  practice.  One  of  the  primary 
goals  of  the  control  strategy  must  be  to  keep  the  tow-line  tension  from  increasing  quickly  as  this  leads  to 
parasail/payload  system  instability.  Furthermore,  wind  speed  sensing  will  be  critical  during  tow-up  and  during  powered 
and  coast  phases  of  flight.  Active  control  during  powered  and  coast  flight,  as  well  as  during  landing,  are  extremely 
important.  Not  only  does  this  facilitate  arriving  to  a  desired  location,  but  also  permits  wind  disturbance  rejection 
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